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We summarize the main motivations to study event shapes at hadron coUiders. In addition we present classes of 
event shapes and show their complementary sensitivities to perturbative and non-perturbative effects, namely jet 
hadronization and underlying event. 



1. Introduction 

Event shapes describe the energy and momentum flow of the final state in high energy collisions. Thanks to the 
synergy of their simplicity and sensitivity to properties of QCD radiation, they are among the most extensively 
studied observables in e+e" and DIS collisions. Specifically, studies of event shapes were used for measurements of 
the strong coupling constant as and of its renormalization group running P|, for cross-checks of the gauge group 
through measurements of the colour factors and, most important, they made it possible to get insights into the 
dynamics of hadronization 0, Q| • 

Despite the great success of these studies in e+e" and DIS, event shapes have been largely neglected at hadron 
colliders, the only (published) exceptions being a measurement of a variant of the broadening by CDF in 1991 [sj and 
of a variant of the thrust by DO in 2002 At hadron colliders, the counterpart of e+e" — > 2 jets or DIS [l+lj-jet 
event shapes are those in high-transverse momentum dijet production. 

Measurements and calculations are more difficult at hadron colliders. Experimentally, one has to deal with the 
omnipresent underlying event and with the impact of a limited detector coverage. Theoretically, the presence of 
four hard QCD-emitting partons at Born level implies complicated patterns of interference between different dipoles, 
which for the first time involve non-diagonal colour structures ^j. Given these difficulties, it is natural to ask oneself 
whether there is anything new that can be learned from event shapes at hadron colliders. In the following I will 
address this question. 

2. New aspects of hadron-hadron event shapes 

It is worthwhile to set at first some terminology and conventions. A given event shape is denoted by V, while v 
denotes its value given a set of momenta, and L = ln(l/w) denotes its logarithm. By convention, event shapes are 
defined in such a way that in the absence of secondary emissions their value is precisely zero. A small value then 
denotes the presence of only soft-collinear emissions. In this region the smallness of the QCD coupling is compensated 
by large logarithmic contributions. Fixed-order perturbative (PT) predictions become unreliable and logarithmic 
enhanced terms have to be resummed to all orders. The state-of-the-art is to resum terms up to next-to-leading 
logarithms (NLL), i.e. a^i", in the exponent of the integrated distribution of w, Ti{v) = ex-p{Lgi{asL) + g2{asL)} |^. 
Resummed predictions are then matched to next-to-leading order (NLO) fixed-order calculations (i.e. O (cks^) relative 
to the Born process), so as to have an accurate description in the whole phase space. 

These accurate perturbative predictions do not account for hadronization effects. Therefore, PT results are sup- 
plemented with non-perturbative (NP) effects describing hadronization. Within the dispersive approach 9j, one 
assumes that, while the perturbative coupling diverges in the infrared (IR) because of the presence of the Landau 
pole, the physical coupling is well defined and integrable in the IR region. One then introduces the average of the 
full QCD coupling in the NP region 




(1) 
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In this framework it is possible to show that leading NP corrections to event-shapes give rise to a shift (5V^^^ of the 
PT distribution which turns out to be universal, i.e. it is given in terms of an observable-dependent, perturbatively 
calculable coefficient, Cv^ times the pure non-perturbative, observable-independent parameter ap: 

With today's techniques, ao cannot be computed from first principles, since this would require the knowledge of 
the QCD dynamics in the IR, this information must therefore be extracted from data. However, the beauty of this 
approach lies in its simplicity and predictive power: having fitted ao from one observable at a given energy, one can 
predict leading l/Q hadronization effects not only at different energies, but also for different event shapes. 

The dispersive approach turned out to be very successful (see 0] and references therein), in some cases even so 
successful that theorists had to revisit their predictions in order to agree with the experimental data [lol |. It is 
however important to keep in mind that the dispersive approach has been tested only with two-jet event shapes, i.e. 
for those observables for which the Born level event is made out of two hard QCD partons and the first non-zero 
contribution starts with configurations with three QCD partons. While there exist theoretical calculations for non- 
perturbative effects in the case of three-jet event shapes in e'^e" or DIS ([2-1-1]) and they predict the same universal 
behaviour for the power corrections [lll | with the same parameter ao, this important prediction has never been tested 
experimentally to date. 

The calculation of power corrections to three-jet observables makes it possible to test the dispersive approach in 
a new regime. Indeed, for two-jet event shapes it is sufficient to assume that NP radiation is distributed uniformly 
in rapidity, as in the Feynman tube model, to recover the predictions of the dispersive approach. For three-jet event 
shapes, this is no longer true, since for radiation in between jets there is simply no natural direction with respect to 
which one can assume rapidity invariance of density of emitted hadrons. Therefore a confirmation, or falsification, 
of such predictions would shed light on the pattern of NP radiation. For a more detailed discussion, see [l^ . 

A different NP approach is based on shape functions |l^, which account not only for leading NP effects, but 
also for sub-leading ones {O {1 / {vQ)^)) . Such an approach is needed in the region of extremely small values of 
V ^ Aqcd/Q- The price paid for the inclusion of higher-order NP effects is the partial loss of predictive power: 
the knowledge of the shape function for one observable is not sufficient to predict hadronization effects for different 
observables. Two important exceptions exist. One is represented by the class of angularities in e+e~ 0|. This 
class of observables is given in terms of a continuous parameter — oo < a < 2 (a = denotes the thrust and a = 1 
the broadening); and for them it has been shown that shape functions are related by a simple scaling rule |l4l llfj]. 
This important theoretical prediction has not yet been tested experimentally. A second class of observables is the 
so-called fractional energy-energy moments (see Appendix 1.2 of [ig)- These observables have the same NP scaling 
rule as the angularities, and, as angularities for a < 1, they are linear in the secondary emissions, which means that 
their resummation is as simple as that of the thrust. 

At hadron colliders, where two hard QCD partons are present already in the initial state, any final-state mea- 
surement will lead beyond the well-tested two-jet regime. Additionally, unlike in e+e^ or DIS, in hadron-hadron 
collisions multijet events are the natural playground for these QCD studies, since they events appear at Born level, 
without any further suppression. Both PT and NP effects are enhanced at hadron colliders, because of the gluonic 
colour charges. For instance in e"'"e~, at 0{as) leading logarithms have an overall colour factor 2C'f- In hadronic 
dijets, for the pure gluonic channel, this is to be compared with 4:Ca- However, standard NP radiation is now always 
accompanied by the presence of the underlying event. While this complicates the studies of the NP corrections, if 
one is able to disentangle the two effects, event shapes can provide useful information about the underlying event. 
Currently, the o nly observable used to extract information about the underlying event is the away-from-jet parti- 
cle/energy flow [l7|, which is however subject to larger theoretical uncertainties than dijet event shapes. In the 
following I will describe how event shapes can be constructed which deliberately enhance or suppress the effect of 
the underlying event, thereby limiting its impact for purely perturbative studies, or enhancing it in order to focus 
on non-perturbative effects. 
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3. Observables 

The way resummed predictions are currently obtained in hadron-hadron collisions is with the automated resum- 
mation tool CAESAR [l^- A current limitation of CAESAR (and of analytical resummations) is that the observable 
should be global 0|, i.e. sensitive to emissions everywhere in phase space. This is because resummations are based 
on the independent emissions approximation S|, which is known to fail to describe all NLL for non-global observ- 
ables, because of configurations with large-angle soft gluons in the unobserved region, which emit large-angle energy 
ordered gluons in the observed region. At hadron colliders this becomes a serious drawback, since detectors can cover 
only a limited rapidity region. Notably no measurements can be performed too close to the beam. This is usually 
parameterized with a maximum accessible rapidity ryo ■ 

Theoretical and experimental requirements thus seem to be in conflict. In spite of this, it was shown in J^J that 
three different classes of observables can be designed to overcome this tension. In the following I will introduce these 
observables and show that they have complementary sensitivities, specifically with respect to the underlying event 
and jet hadronization effects. 



T„,, ^ ; . (4) 



3.1. Directly global observables 

To define directly global observables one selects events with two large transverse momentum jets and defines 
observables similar to e+e^, but purely in the transverse plane. For instance the directly global transverse thrust is 
given by 

Tj_^g = max ' 1^ I — , rj_,g = 1 - Tj_,g , (3) 

where p±^i denotes the transverse momentum with respect to the beam of parton i and the directly global thrust 
minor is defined as 

Similarly to event shapes, one can define three-jet resolution variables. In the longitudinal invariant exclusive kt 
algorithm one defines the distance between particle i and the beam B and between two particles i,j 

d,B = q\i , d,j = min{(?i,, q\j} ((77^ - if^) + (0, - (ji^f) (5) 

and then recombines the pair with smallest distance. Finally, one defines 

where d*^"-' is the smallest distance when n pseudo-particles are left over after recombination and E±i^2 denote the 
transverse energies of the two leading jets. As for dijet event shapes, 2/23 1 denotes a two-jet-like configuration 
(see Fig.[T}. 

For directly global observables, measurements should be carried out as forward as possible (up to rjo ^ 3.5 at the 
Tevatron and 770 ~ 5 at the LHC). It has been shown |2(| that measurements in the unobserved region are negligible 
at NLL as long as the observable is not too small. More precisely, given an observable which, in the presence of a 
single soft emission collinear to leg £, behaves as 

^-rf£(|)"e"'"V(0), (7) 

the formal limit for NLL predictions to be valid reads 

\ii{X-V)>-{a + b)-rjo, 6= min{fei,62}- (8) 

Here X denotes a rescaling of the argument of the logarithm chosen so as to cancel the average value of degi{(j)) in 
the resummed exponent g2{asL), see 0,0] and 1,2 label the incoming legs. The strategy is then to neglect tiq in 
theoretical predictions and to check a posteriori which portion of the tail is beyond control. 
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J/23 ~ 1 



2/23 < 1 





Figure 1: Illustration of a three-jet event, i/23 ~ 1 (left) and of a two-jet-like event, 7/23 ^ 1 (right) 

3.2. Recoil enhanced observables 

The definition of recoil enhanced observables proceeds as follows: define a central region C, e.g. \ri\ < rji-nnK ~ O (1), 
which should contain the two jets, see Fig. [3 Then define central observables using particles only in the central 
region. For instance, the central transverse thrust is given by 



Tl,c = max =^ — — 



l-T 



_L.C ■ 



(9) 



Similarly one defines a central thrust minor, or three-jet resolution variable, as in eqs. Q and ^ but restricting the 
measurement to the central region. 

Transverse-momentum conservation ensures that 



1 



i£C 



1 



(10) 



Central observables can then be made global by adding the recoil term 7?,j_^c (or a power of it), e.g. the recoil enhanced 
transverse thrust, the thrust minor and the three-jet resolution parameter read 



m.C 



(11) 



Note that to ensure continuous globalness |21| (roughly speaking that the transverse-momentum dependence be 
independent of the emission's direction) one has to add a second power of the recoil term to 1125,0 ■ 

Observables in this class are global despite the fact that measurements are truly restricted to a central region. 
This is because they have an indirect sensitivity to forward emissions through recoil effects. This mechanism was 
already well known from some DIS observables 



;y_t( 





C 

Figure 2: The central region C, containing the jets, and the forward rapidity cut r]o. 
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3.3. Observables with exponentially suppressed forward terms 

The last class of observables discussed is the one with exponentially suppressed forward terms. Also in this case one 
starts defining central observables, as in eq. One then introduces the mean transverse-energy weighted rapidity 
rjc of the central region 



■nc 



and defines an exponentially suppressed forward term as follows 

^c^TT^E?^^^""""'"'' (13) 

where the total transverse momentum Q±.f, of the central region is defined in eq. H10|l . With these elements, one can 
make the usual central observables global by adding the exponentially suppressed forward term £q (or a power of 
it to ensure continuous globalness). For instance, the exponentially suppressed transverse thrust, thrust minor and 
three-jet resolution parameter are given by 

= Tj_,c + ' Tm,£ = T„t.c + £c , y23,e = 2/23,c + £c ■ (14) 

Additionally, one can consider observables that are more suitably defined only in a restricted (central) region. One 
first separates the central region C into an up part Cu, consisting of all particles in C with p± ■ nT,c > 0, and a down 
part C^), consisting of all particles in C with p± ■ ut.c < 0. 

One then defines, in analogy to e"'"e~ |i22i |. the normalized squared invariant masses of the two regions: 



Pxfi^j^(Y.1^ ' X^U.D, (15) 

from which one can obtain a (non- global) central sum of masses and a (non-global) heavy-mass: 

Psfi = Pufi + PDfi , PH,c = max{/9c/,c,Pr>,c} , (16) 

together with versions that include the addition of the exponentially suppressed forward term: 

PS£ = Psfi + £c > Ph,£ = PHfi +£c- (1"^) 

This separation into up and down regions can also be used to define boost-invariant jet broadenings. One first 
introduces rapidities and azimuthal angles of axes for the up and down regions: 

Cx V TT r\ n o^ 

VX.C = , 0X,C = -^^ , X = U,V , (iSJ 

and defines broadenings for the two regions: 

Bx,c = I^^^Ji^^ - ^x,c? + (0» - 0x,c)2 , X = U,D. (19) 

The central total and wide-jet broadenings are then given by 

Brfi = Bufi + Bofi , Bwfi = max{Bufi, Bofi} ■ (20) 

Exponentially suppressed, global observables are then obtained by adding the forward term £q 

Bt,£ = Brfi + J Bw,£ = Bwfi + ■ (21) 
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V 


ae (inc) 


be (inc) 


ai (out) 


be (out) 


T-L,9 


1 





1 


1 


T 


1 





1 





1/23,9 


2 





2 






V 


ae (inc) 


be (inc) 


ae (out) 


be (out) 


PH,E 


1 


1 


1 


1 


Bw,e 


1 


1 


1 






Table I: Coefficients a and be for some global (left) and exponentially suppressed (right) observables as they appear in eq. lO, 
given for incoming and outgoing legs. 

4. Sample NLL distributions from caesar 

To illustrate different features of resummed predictions, we present here some sample distributions as obtained 
with CAESAR. We select dijet events at the Tevatron Run II regime, -y/s = 1.96 TeV, by running the longitudinally 
invariant kt algorithm and by requiring the presence of two jets in the central region (I77I < 0.7), the transverse 
momentum of the hardest jet satisfying E± > Ej__inin = 50 GeV (unless differently specified). We define the central 
region C \ri\ < ?7inax = l-l- This ensures that the two jets are well contained in C. 

Most of the properties of distributions can be understood in terms of the coefficients a and bi appearing in the 
parametrization of the observable given in eq. Q. This is because the leading Sudakov effect is given by 

E(l/)^e5#G..i^+-, Gi2 = --E4V' (22) 

a ^ — ' a + oi 

i 

so that it is determined only by the colour charges Cg {Cg = Cp, Ca if parton £ is a quark/gluon) and by the values 
of a and bf . The point up to which the resummation is under NLL control, given in eq. ©, also depends on the 
values of a and bg (although it also depends on the values of digi{ip) through X). 

In TableQlwe report the values of the coefficients a and bi for various observables as computed by CAESAR (despite 
the presence of four legs it is sufficient to distinguish only between in/outgoing legs). 

In Fig. Owe show pure NLL resummed differential distributions for various global observables, the thrust, thrust 
minor, and three-jet resolution parameter, defined in eqs. Q, and © for various partonic channels. The figure 
shows that the thrust minor is peaked at higher values of u, while the peaks of y23 lie at lowest values. This is 
expected since for these observables the position of the peak depends only on the coefficients of the total colour 
charge of in/outgoing partons, Ci„/out (specifically here G12 = -(2Cin + Cout), -2(Cin + Cout),-(Cin + Cont) for 
the thrust, thrust minor, and three-jet resolution, respectively). The figure also shows that tails of quark-channel 
distributions generally extend to lower values of v compared to channels dominated by gluons. This can again be 
inferred from the value of G12 for the given channel; it simply refiects the fact that quarks radiate less then gluons. 
The region to the left of the thick line is beyond control of the NLL resummed predictions, given in eq. l|HJl. Note 




Figure 3: NLL resummed differential distributions for the directly global thrust (left), thrust minor (central) and three-jet 
resolution parameter (right), for various partonic channels. 
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ln(V) ln(T^_R) ln(T,^_p,) 

Figure 4: Total differential distribution for exponentially suppressed wide broadening and heavy jet-mass for two valued of 
£'±,min (left). Differential distributions for the recoil enhanced thrust (central) and thrust minor (right) for various partonic 
channels. 



that X h for the thrust, while X = 1 for the thrust minor and 1/23, so that for these last two observables the cut 
is located simply at 770 and 2r]Q respectively. We see that for the thrust a small part of the tail extends beyond the 
cut, while for the other two observables almost all the distributions arc under pcrturbative control. 

Finally, one should remark that at higher values of v NLL distributions become negative. This unphysical behaviour 
will be corrected by the fixed-order predictions once the matching with fixed order results is carried out. 

In Fig. 21 (left) we show total distributions, summed over all partonic channels, for two exponentially suppressed 
observables, a variant of the heavy jet-mass, and the wide broadening defined in eqs. H17|l and (|21(l for two values of 
-E_L,min- The two observablcs have a different leading Sudakov suppression (6*12 = — (Cin + 2Cout), — (Cin + Cout) for 
the broadening and the jet mass respectively). Compared to directly global observables, one can see a remarkable 
improvement as far as the position of the cut is concerned, due to the exponentially suppressed term [bp = 1 for 
incoming legs). 

One can also see the effect of changing the transverse energy cut: increasing the cut, the distribution is peaked at 
lower values. This is due to two different mechanisms: with higher transverse energy cuts one probes the coupling 
at higher scales and the parton densities at higher x, where quark channels become more important. 

Finally, in Fig. 0] (centre and right) we show the distribution for the recoil-enhanced transverse thrust and thrust 
minor defined in eq. Since these observable measure only very central emissions, the 770 cut has no effect. 

However, the distributions have a different cut since the resummation breaks down at some value of v (which 
generally depends on the channels too) . This breakdown is well understood [3 , and it is a characteristic of all 
those observables for which contributions from multiple emissions can cancel one another. What happens is that, for 
sufficiently small v, it becomes more likely to suppress the value of v via cancellations between emissions rather than 
by Sudakov suppression. This leading logarithmic (LL) change of behaviour cannot be accounted for by the NLL 
resummation. Points where the NLL resummation is not reliable have been excluded from the plots. The effect of 
this can be seen in the cut of the tails of the distributions. 



5. Complementarity between the observables 

The above observables have complementary properties, sensitivities and advantages/disadvantages, which makes 
them a very rich field of activity. To be more precise, directly global observables are simple, being defined analogously 
to e^e~ observables, just purely in the transverse plane. However, the impact of 770 is not always negligible and this 
can be checked only a posteriori. Observables with recoil term measure only central emissions. The recoil term can 
however be subject to large cancellations, making an accurate determination of it problematic. Additionally, the 
Sudakov effect competes with a cancellation in the vectorial sum in the recoil term, causing the resummation to break 
down at some point. Observables with an exponentially suppressed forward term have a very low sensitivity to the 
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Event shape 


Impact of i]o 


Resummation 
breakdown 


Underlying event 


Jet hadronization 




tolerable* 


none 


~ Vo/Q 


~1/Q 


Tm,g 


tolerable 


none 


~ Vo/Q 






tolerable 


none 






T±,e, px.e 


negligible 


none 


~1/Q 


~1/Q 


Bx,e 


negligible 


none 


~1/Q 




Tm,E 


negligible 


serious 


~1/Q 




V2-i,e 


negligible 


none 


~1/Q 




T^,TZ, PX.TZ 


none 


serious 


~1/Q 


~1/Q 


Tm,TZ, Bx.n 


none 


tolerable 


~1/Q 


~ 1/(VS7Q) 


y23,TZ 


none 


intermediate* 







Table 11: Summary of the complementary properties of various observables. Entries marked with a * are subject to uncertainty 
at present. 



forward region; they are thus the most promising compromise. Yet experimental measurements are complicated by 
the fact that one needs to merge measurements from different detectors (tracking information in the central region, 
where a good resolution is needed, calorimeter information in the forward region, where an "average" information 
on the momentum flow is sufficient). 

These theoretical complementarities are summarized in the first two columns of Table IIII for variants of the thrust 
(tj^), thrust-minor (Tm), broadenings (Bx)-, jet-masses (px), and three-jet resolution parameter (y23)- Additionally 
the table shows the sensitivity of these observables to the underlying event and the hadronization of the jets. In the 
following we will explain the listed properties. 

Let us begin by considering effects of jet hadronization. Studies of e^e~ observables showed that, for those that 
are exponentially suppressed in rapidity, the power corrections amount to a rigid shift cx l/Q ^2^, while those that 
are uniform in rapidity power corrections in addition lead to a squeeze of the distribution which is cx Xji^^falQ) [lol |: 

V ^ ^ ^ - cy(w)ao 



The origin of the observed difference is the following. For an observable V ^ which behaves as ^ ~ ^fviv)i for a 
fixed perturbative configuration, the contribution to the power correction from parton with colour charge can be 
written schematically as 

AT/ f *-Lj CiOsikjJkjL , . , , 

^^-^ / -[—'^^ -rrlviri), (23) 

where the integration is over the region where the gluon is non-perturbative. 

For observables which are exponentially damped in rapidity, the rapidity integration can be extended up to in- 
finity. The rapidity and transverse-momentum dependence then completely factorize: the rapidity part goes into 
the observable dependence coefficient cy = J dr]fv{ri), while the transverse- momentum integral gives rise to ao, the 
average of the coupling in the IR region, below some merging scale given in eq. One obtains 

SV cx cvao ^ . (24) 

For observables which are uniform in rapidity, the rapidity integral cannot be extended up to infinity and PT 
radiation provides an effective cut-off in the rapidity integration. This cut-off is generally given by the angle between 
the direction of hard emissions and the reference directions used for the measurements (e.g. thrust axis in 2-jet events, 
event-plane for near-to-planar 3-jet events). This angle depends logarithmically on the transverse momentum of the 
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hard emitting parton qt. Averaging this contribution over the Sudakov PT recoil distribution e~"^'" iit/Q) amounts 
to a singular correction cx 

This same pattern is expected to hold in hadron-hadron collisions as well. The parametrization of the observable 
given in eq. lO, which is established automatically by CAESAR, allows one to determine the rapidity behaviour of 
the observable (i.e. the value of the bg) and therefore the type of power correction. Specifically, observables with 
bg — 1 for outgoing legs will have power corrections 1/Q, while those with bg — for outgoing legs will have power 
corrections l/{^/a^Q). Finally, one sees from the table that power corrections for the j/23 jet-resolution parameters 
are estimated to be cx ^yxilQ- The reason is that, contrary to the other observables, which are linear in kt^ 2/23 is 
quadratic, therefore the effect of the power correction can be estimated via 

<52/23 = y23,NP - J/23,PT oc J — - — dr]as{kj_)\ ^ — I (25) 

f^i jlNP i.PTl,NP 

-Q—sik^)-Q^ - ■ (26) 

Finally, let us motivate how the underlying event affects measurements of these observables. We assume that 
particles from the underlying event are distributed uniformly in rapidity |27| : therefore the effect of the underlying 
event to an observable V, parametrized according to eq. Q, can be estimated by 

(^(U.E.)^^ - /" J dk^^^^d,^^ ['^y e~^^^^Mg,M) - ^-^^^idM) /" '^^e""-!"! . (27) 

If we consider observables for which bg — for incoming legs, the underlying event will be enhanced by rjQ, while for 
observables with bg ~ 1, the rapidity integral can be safely extended up to infinity and the result is 770-independent. 
For observables quadratic in k±, like j/23, taking into account the interplay between PT and NP emissions amounts 
to an effect that is observable-dependent and goes again as y^yis/Q. Finally, there are observables for which there 
is a suppression coming from the vanishing of {gi^2{<t>))- This suppression is the reason why the contribution of the 
underlying event does not have an t^q enhancement in the case of Tm,£, despite the fact that bi/2 = 0. The different 
ways in which jet hadronization and underlying event affect the event shape distributions should make it possible to 
disentangle the two effects and study them separately. 



6. Conclusions 

We showed that different observables in hadron-hadron collisions have complementary sensitivities and properties. 
This makes them a powerful tool to investigate properties of QCD radiation, in particular the effects due to jet 
hadronization and to the underlying event. These studies are feasible thanks to the automated resummation tool 
CAESAR. However, before phenomenological studies can be carried out, resummed predictions need to be matched 
with NLO results. In hadronic collisions a channel-by-channel matching requires the flavour information [2^ for the 
fixed-order predictions, which is unfortunately not currently available to an external user of hadron-hadron NLO 
programs [291 l3C||. Work in this direction is in progress. 
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